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Abstract 

We update indirect constraints on Electro-Weak Dark Matter (EWDM) considering the 
Sommerfeld-Ramsauer-Townsend (SRT) effect for its annihilations into a pair of standard model 
gauge bosons assnming that EWDM accounts for the observed dark matter (DM) relic density for 
a given DM mass and mass gaps among the multiplet components. For the radiative or smaller 
mass splitting, the hypercharged triplet and higher multiplet EWDMs are ruled out up to the 
DM mass ~ 10 - 20 TeV by the combination of the most recent data from AMS-02 (antiproton), 
Fermi-LAT (gamma-ray), and HESS (gamma-line). The Majorana triplet (wino-like) EWDM can 
evade all the indirect constraints only around Ramsauer-Townsend dips which can occur for a tiny 
mass splitting of order 10 MeV or less. In the case of the doublet (Higgsino-like) EWDM, a wide 
range of its mass > 500 GeV is allowed except Sommerfeld peak regions. Snch a stringent limit on 
the triplet DM can be evaded by employing a larger mass gap of the order of 10 GeV which allows 
its mass larger than about 1 TeV. However, the future CTA experiment will be able to cover most 
of the unconstrained parameter space. 
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I. INTRODUCTION 


The most simplistic candidate for DM would be a neutral component of an SU(2)^ x U(l)y 
multiplet added to the Standard Model (SM) |I1[2], which is dubbed as EWDM. In the case of 
a fermionic EWDM, its physical properties are completely determined by the gauge charge, 
dark matter mass, and mass differences between the multiplet components. As is well-known, 
the non-perturbative effect plays a crucial role when the DM is slowly moving and the DM 
mass is larger than the force-carrying boson mass [IlEHZ]- Furthermore, this effect includes 
not only the usual Sommerfeld enhancement but also the Ramsauer-Townsend suppression 
which are more apparent for larger DM mass or smaller mass gaps [7]. 

In this paper, we revisit the non-perturbative effect of EWDM which is strongly con¬ 
strained by the recent indirect detection data on anti-proton flux by AMS-02 [S], gamma-ray 
measurement from Milky Way (MW) satellite dwarf galaxies by Fermi Large Area Tele¬ 
scope (Fermi-LAT) [9], and gamma-ray line searches by Fermi-LAT [10] and High Energy 
Spectroscopic System (HESS) [H]. The first two measurements put bounds on the leading- 
order annihilation processes: —)■ WWjZZ^ and the last two on the loop-induced ones: 

—)■ I^Z. Remark that the updated Fermi-LAT search on gamma-ray lines covers the 
dark matter mass up to around 500 GeV complimenting the previous HESS search range of 
500 GeV - 25 TeV. Related studies have been made for the case of the wino or wino-Higgsino 
dark matter in a supersymmetric theory [l 2 Hn|. 

The annihilation of two neutral particles into 'y + X is a radiative process and may 
be subject to a large correction due to the resummation of electroweak Sudakov logarithms 
02 log^(’^DM/'^w)- Refs. [miTn], the authors studied the Sudakov resummation effect for 
the annihilation rate of the exclusive process x^X^ 111 ^ind found that the exclusive 
rate is reduced by up to a factor of 2 - 3 compared to the tree-level plus Sommerfeld 
enhancement calculation. On the other hand. Refs. [ 20 ] treated systematically the semi- 
inclusive annihilation rate into the hnal state 7 -|- X within the resolution of the experiment 
to hnd that the effect of higher order correction is very limited as the semi-inclusive rate 
changes by only 0(1%) at tudm = 3 TeV. Thus, we will simply use the leading-order 
annihilation cross section of EWDM including only the SRT effect. 

For our study, we will consider three specific examples of fermionic EWDM in the lowest- 
lying multiplets: a vector-like (Dirac) doublet with a hypercharge Y = ±1/2 (Higgsino- 
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like), a (Majorana) triplet with Y = 0 (wino-like), and a vector-like (Dirac) triplet with 
Y = ±1. Note that a certain symmetry like Z 2 has to be imposed for the stability of 
these EWDM candidates. They are assumed to form 100 % of the observed DM density in 
wide ranges of the DM mass and mass gaps among the multiplet components. As shown 
in Ref. [7], the feature of Sommerfeld peaks and Ramsauer-Townsend dips in the non- 
perturbative annihilation cross section depends sensitively on those mass parameters, which 
leads to an interesting impact on the indirect detection constraints. 

This paper is organized as follows. In Section II, we review the non-perturbative effect, 
the SRT effect for annihilations of EWDM summarizing the results in Ref. [7j. In Section 
III, we discuss the existing indirect constrains on annihilations of EWDM into SM gauge 
bosons based on various cosmic-ray measurements, then conclude in Section IV. Finally, in 
Appendix A, we present the explicit forms of the potential and scattering matrices described 
in Section II. 


II. EWDM AND SOMMERFELD-RAMSAUER-TOWNSEND EFFECT 

Let us hrst make a brief summary of the non-perturbative effect on the EWDM annihi¬ 
lation [7]. The doublet EWDM (Higgsino-like) consists of a Dirac fermion pair denoted by 
Xu = iXu) Xu) Xd = iXdy Xd) chiral representation. The lighter linear combination 

of xt Xd candidate and the mass splitting can come from (effective) non- 

renormalizable operators after the electroweak symmetry breaking. The wino-like EWDM 
multiplet, a triplet with V = 0 is denoted by y = (x"*", x°, x~) having only one Majorana 
neutral component. Finally, the triplet EWDM multiplet with Y = ±1 consists of two chiral 
fermions, Xu = ixt^,xt:Xl) and Xd = (x°,X7Wd")- Recall that the electroweak one-loop 
correction generates a mass splitting of order 0.1 GeV between the multiplet components 
[2]. Together with the above-mentioned (tree-level) contribution, this one-loop correction 
can make arbitrary mass gaps as assumed in our analysis. 

The non-perturbative effect in the EWDM annihilation arises from the exchange of the 
electroweak gauge bosons which mixes together the two-body states of the multiplet com¬ 
ponents. In the case of the doublet EWDM, there are three states formed by the charged 
(Dirac) component and two neutral (Majorana) components: xiXdy XiXi: ^md where 

X^ denotes the dark matter component. For the wino-like EWDM, there are two two-body 
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states: x~^X and The triplet EWDM with y = ±1 has four two-body states: Xu^Xd > 
XuXd^ and 

For the calculation of the non-perturbative effect, we apply the analysis of Refs. HI Eli 
in which the Green’s functions gij corresponding to the transition between the two-body 
states i and j are shown to follow the Schrodinger equation: 

- J ^ + Vik{r)gkj{r) = Kgij{r) , ( 1 ) 

^DM cir 

with the boundary condition 5 'p( 0 ) = 6ij and dgij{oo)/dr = — Vii{oo))gij{oo). 

Here K = is the total kinetic energy of the two initial dark matter particles in the 

annihilation process, where (3 is the DM velocity. Then, the non-perturbative annihilation 
cross section of the dark matter is 

avix\^ ^ AB) = 2dodljTf , (2) 

where = gQj{oo) and v = 2/3 is the relative velocity between the two incident DM 
particles. Here A and B run over the gauge bosons (hF+, IT", Z, 7 ), and the gauge boson 
final states AB can be W~^W~, ZZ, jZ, or 77 . The explicit forms of the potential matrix 
Vij in Eq. (j^ and the scattering matrix Td^ in Eq. (|^ are collected in Appendix [A| 

The non-perturbative calculation of the EWDM annihilation exhibits not only the usual 
Sommerfeld enhancement with resonance peaks, but also a vanishing cross section realizing 
the Ramsauer-Townsend effect for particular choices of the model parameters. These effects 
are particularly sensitive to the the mass splittings between the dark matter and the charged 
components of the multiplet, and can appear even for the DM mass below the TeV scale 
when the mass splitting is reduced to (9(10) MeV or less. One can find a detailed analysis 
on non-perturbative effects on the annihilation cross section of the EWDM in Ref. [7j. The 
appearance of the Ramsauer-Townsend dips is of a particular interest as it can allow the 
EWDM to evade strong constraints from various indirect detections. 

III. CONSTRAINTS ON EWDM FROM INDIRECT DM SEARCHES 

A. Annihilation into WW & ZZ final states 

We first present constraints on EWDM from indirect signals by its annihilation into 
WW/ZZ final states. Annihilations of EWDM are expected to produce W/Z bosons plen- 
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tifully since EWDM is charged under the SU(2) symmetry of the SM. Fragmentation of the 
produced WjZ bosons leads to sizable contributions to the antiproton flux and the contin¬ 
uum photon spectrum which are detectable in cosmic-ray measurement experiments such as 
AMS-02 and Fermi-LAT. 

1. Constraints from AMS-02 antiproton flux measurements 

Antiproton production from DM annihilations into WW and ZZ channels is constrained 
by the precise measurements on antiproton flux by AMS-02 [S] . The antiproton contribution 
from DM annihilation should be summed to the secondary antiprotons, produced by colli¬ 
sions of energetic cosmic-rays with the interstellar medium, which account for the bulk of the 
observed flux. Ref. |2T] provided constraints on DM annihilation cross sections for various 
hnal states including WW based on the recent antiproton to proton ratio measurements by 
AMS-02 [H]. The authors of Ref. [21] assumed the Einasto DM halo prohle and the MED 
propagation parameter set proposed in Ref. [22|- They also showed that for /udm > 100 
GeV the limits are ~ 2 times weaker for the Burkert prohle and 2 — 3 times more stringent 
for the MAX propagation modelj^ In our study, we use the 2a exclusion bound on the WW 
channel for the Einasto prohle with the MED propagation model in Ref. [2l] as the bound 
on the total cross section av'^'^ -|- av^^ = av^o^o^^y+iY- -\- av^o^o^zz since the antiproton 
yields per annihilation from two hnal states WW and ZZ are almost undistinguishable. The 
AMS-02 antiproton hux limit on av^^ -\-av^^ is shown as a red-thin dotted line in Figures 

mm 

2. Fermi-LAT continuum photon constraints: Milky Way satellite dwarf galaxies 

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most promising 
targets for the DM indirect detection via gamma-rays since they are highly DM domi¬ 
nated objects with a relatively short distance from the Earth. In EWDM annihilation, 
the continuum photons originate mostly from fragmentation of hadronic hnal states in the 
WWjZZ processes which are strongly constrained by the gamma-ray measure- 
^ In a number of recent papers based on Galactic synchrotron emission |2dH2B) and cosmic-ray positrons |241 
HZlllH], it has been pointed out that the thin halo model in the MIN propagation scheme is seriously 
disfavored. 
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ments from the MW satellite dwarf galaxies [9]. The total mass of a dwarf galaxy within 
the half-light radius and the integrated J-factor, the line-of-sight integral through the DM 
distribution integrated over the solid angle for a region of interest, have been found to be 
pretty insensitive to the change of a DM density prohle [29H3T] . Thus in Ref. [9] , the bound 
from satellite dwarf galaxies is obtained based on 6 years of Fermi-LAT data processed with 
the Pass 8 event-level analysis just assuming that the DM distribution in dwarf galaxies 
follows a Navarro-Frenk-White (NFW) prohle |32]. In our analysis, we take the Fermi-LAT 
continuum photon limit on {+av^^) at 2(j level for the MW satellite dwarf galaxies 

from Ref. [9] which is shown as a red-thin dashed line in Figures 


B. Annihilation into 77 &: 7 Z final states 

In this subsection, we provide indirect constraints on EWDM by its annihilation into 
77 / 7 Z final states. Monochromatic photons arise from the one-loop processes 
77 / 7 Z. Such a line signature would be quite easily distinguished from astrophysical photon 
sources since in most cases they produce continuous spectra. Now, the annihilation cross 
sections for the processes 11 ^ire already in tension with Fermi-LAT and HESS 

searches for line-like spectral features in the photon spectrum. 


1. Constraints from Fermi-LAT photon line searches 

Very recently, the Fermi-LAT collaboration has reported a constraint on a DM annihila¬ 
tion cross section based on updated searches for spectral line signatures in the energy range 
200 MeV - 500 GeV from around the Galactic Genter (GG) using 5.8 years of data repro¬ 
cessed with the Pass 8 event-level analysis [10] . In the analysis, they searched spectral lines 
expected from dark matter annihilation for four signal regions of interest (ROIs), selected to 
optimize sensitivity to different dark matter halo profiles: NFW with 7 = 1 or 1.3, Einasto, 
and isothermal. In Ref. Hi, it has been shown that the bound for tudm ^ 100 GeV is just 
2-4 times weaker even for the isothermal profile compared to the bound for the Einasto 
profile since the Fermi-LAT has measured gamma-rays from all the sky and thus can find 
the corresponding optimized ROI for each DM profile. In this study, we use the 2a upper 
limit on av"^^ for the Einasto profile with the local DM density po = 0.4 GeV/cm^ as the 
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bound on the total cross section + ^crv'^^ = av^o^o^^j + weighted by the 

number of photons for each hnal state. In Figures[T]-[^ we plot the constraint for the Einasto 
prohle on crn'’"^ + as a blue-thick dotted line. 

S. Constraints from HESS photon line searches 

In Ref. m, upper limits on line-like gamma-ray signatures in the energy range 500 
GeV - 25 TeV are provided using the data collected by the HESS, which complement 
the limits obtained by the Fermi-LAT at lower energies [TOj. In the analysis, the HESS 
collaboration assumed the Einasto DM halo prohle with po = 0.4 GeV/cm^. However, the 
HESS collaboration has searched gamma-ray lines for only one ROI, a 1° radius circle around 
the GG, compared to the Fermi-LAT. Thus, the HESS limit for the Einasto prohle can be 
weakened by about two orders of magnitude for a more cored prohle, the isothermal prohle. 
In our analysis, we use the 2 ct limit from the HESS spectral line search for the Einasto prohle 
as a representative constraint on av'^'^ + which is depicted as a blue-thick dashed line 

in Figures [T]-|^ 

C. Summary of indirect constraints on EWDM 

We plot DM annihilation cross sections and + av^^ as blue-thick and 

red-thin solid lines for the doublet, Majorana triplet, and hypercharged triplet {Y = ±1) 
EWDMs in Figures and respectively. For each type of EWDM, we calculate av'^^ 
+ \crv'^^ and as a function of DM mass tudm including the non-perturbative 

ehect for two representative values of the typical mass splitting of 0(0.1) GeV arising 
from the electroweak one-loop correction (top-panel) and 0(10) MeV (bottom-panel). For 
larger mass gap, the SRT ehect becomes much weaker for a given DM mass. In order to 
see this ehect, we additionally show the results for the V = ±1 triplet EWDM with the 
mass gap of 10 GeV in Figure]^ All the relevant indirect constraints on and 

av^w -|- (jv^^ are shown by blue-thick and red-thin curves, respectively. 

• The limit from AMS-02 antiproton flux measurements: We plot the 2a limit 

on av^'^ + av^^ from AMS-02 antiproton flux measurements as a red-thin dotted 
line. 
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Higgsino-like, 5m+=34I MeV, <5mjv=0.2 MeV 
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FIG. 1. Constraints on -l- ^v^o^o^zz (red-thin curves with the scale on the right) 

and -|- ^av^o-^o^^z (blue-thick curves with the scale on the left) for the Higgsino-like 

EWDM with 6m+ = 341 MeV (top) and 8 MeV (bottom) for 6mN = 0.2 MeV. In each panel, the 
solid line is the calculated DM annihilation cross section including the non-perturbative effect. The 
red-thin dotted and dashed lines show the upper limits obtained from the AMS-02 anti-proton flux 
data analysis m and Fermi-LAT gamma-ray measurements from MW dwarf spheroidal satellite 
galaxies [9] . The blue-thick dotted and dashed lines are constraints from line-like photon signature 
searches by Fermi-LAT [TO] and HESS [TT] . respectively. The CTA sensitivities on the WWjZZ |33| 
and 7 A [341135] channels are presented as red-thin and blue-thick dot-dashed curves, respectively. 
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Wino-like, 6m+=l66 MeV 
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FIG. 2. Constraints for the wino-like EWDM with = 166 MeV (top) and 6 MeV (bottom). 
Each line is the same as Eigur^ 


• The limit from Fermi-LAT continuum photon searches: For av^'^ + av^^, the 
upper region of the red-thin dashed curve is excluded by the limit from the Fermi-LAT 
continuum gamma-ray searches for the Milky Way satellite dwarf galaxies at the 2a 
level. 
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Triplet Y=l, 5m++=1.4 GeV, ^m+=525 MeV, <5mft.=0.2 MeV 




FIG. 3. Constraints for the triplet EWDM with Y = ±1. 6m+ = 525 MeV (top-panel) and 15 
MeV (bottom-panel) for = 1.4 GeV and 6m]y = 0.2 MeV. Each line is the same as Eigur^ 

• The limit from Fermi-LAT photon line searches: For , the 2a 

exclusion limit from Fermi-LAT spectral line searches for around the GC is shown as 
a blue-thick dotted curve. 


• The limit from HESS photon line searches: The HESS gamma-ray line signature 
search limit on av'^'^ -f ^av'^^ is plotted as a blue-thick dashed curve. 
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“DM [GeV] 


FIG. 4. Constraints for the triplet EWDM with Y = ±1 for = 20 GeV, (5m+ = 10 GeV, 

and 5mi\f = 0.2 MeV. Each line is the same as Figur^ 

• The future sensitivity of CTA; In near future, the remaining parameter regions 
for each EWDM will be probed by various upcoming cosmic-ray observations such as 
CTA |36] and GAMMA-400 [37] . The CTA sensitivity with a 500 h time exposure on 
(jyWW _|_ (jyZZ jggj jg plotted as a red-thin dot-dashed line. For 5 h of GC observation 
with CTA, the upper limit on av'^'^ -|- [Ml |35] is depicted as a blue-thick dot- 

dashed curve. 

As stated in Introduction, each EWDM is assumed to account for 100 % of the observed 
DM relic abundance in wide ranges of the DM mass and mass gaps. In the case of the doublet 
(Higgsino-like) EWDM with a radiative mass splitting 6m+ = 341 MeV, the DM mass larger 
than about 500 GeV is allowed except a narrow Sommerfeld peak region at 7 TeV. For a 
smaller mass gap = 8 MeV, the hrst Sommerfeld peak moves down to about 1 TeV and 
more peaks appear at lower DM masses. These peak regions are excluded by either the Fermi- 
LAT gamma-ray data from the dwarf galaxies or HESS gamma-line data. The Majorana 
triplet (wino-like) EWDM is stringently constrained for the whole range of masses up to 
~ 10 TeV for the typical mass splitting of 166 MeV by the AMS-02, Fermi-LAT, and HESS 
results. However, the indirect constraints on the wino-like EWDM could be evaded around 
Ramsauer-Townsend dips which can occur for a very small mass splitting, 0(10) MeV or 


11 







less. For the case of the triplet EWDM with Y = iblj^the indirect constraints become much 
more stringent to rule out for masses up to ~ 20 TeV by a combination of AMS-02, Fermi- 
LAT, and HESS limits. It is also difficult to dodge the indirect constraints by arranging 
Ramsauer-Townsend dips with smaller mass splittings as, contrary to the Sommerfeld peaks, 
the Ramsauer-Townsend dips for + av^^ and av"^"^ + ^crv'^^ do not coincide with each 

other as shown in the lower panel of Figure 3. This behavior is more apparent in the case of 
the quintuplet EWDM which exhibits non-overlapping narrow dips even with the radiative 
mass splitting |3], and thus is completely ruled out by the combination of all the indirect 
constraints. As shown in Figure taking a large mass splitting of 0(10) GeV, the triplet 
EWDM with Y = ±1 can be safe from all the current indirect constraints in the mass range 
of rriBM ^ a couple of TeV except around the Sommerfeld peaks. Similarly, the Wino-like 
EWDM can escape from the current indirect limits for ttidm ^ 1 TeV. However, CTA [36] 
will be able to probe most of remaining parameter regions even for the Higgsino-like EWDM 
and/or a large mass splitting of 0(10) GeV. 

IV. CONCLUSION 

In this paper, we have investigated indirect constraints on EWDM considering the SRT 
effect for its annihilations into a pair of SM gauge bosons which are sensitive to the size of 
mass splitting among the multiplet components. Assuming that EWDM accounts for the 
observed DM relic density for a given DM mass and the radiative (or smaller) mass splitting, 
we found that the triplet with Y = ±1 (and higher multiplet) EWDM is ruled out up to 
moM ~ 10 - 20 TeV by the combination of current limits from AMS-02 (antiproton), Fermi- 
LAT (gamma-ray), and HESS (gamma-ray line) measurements, disregarding a potentially 
strong DM halo prohle dependence of the HESS limit. In the case of the Majorana triplet 
(wino-like) EWDM, there is a chance to dodge the indirect constraints around Ramsauer- 
Townsend dips only with a tiny mass splitting 5m_|_ < 0(10) MeV. On the other hand, 
the Higgsino-like EWDM is excluded just for DM masses less than ~ 500 GeV and around 
Sommerfeld resonance peaks. Such a stringent limit can be weakened signihcantly if a large 
mass splitting of 0(10) GeV is employed. The indirect constraints could be evaded even for 
^ It is advocated as the unique candidate for asymmetric EWDM in Ref. jSSj, which is however stringently 
constrained by DM indirect detection limits as shown in this study. 
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the wino-like and the F = ±1 triplet EWDM in the mass range of itidm > Cl(TeV) except 
aronnd the Sommerfeld peaks which occnr at larger DM masses. However, the nnconstrained 
parameter regions will be mostly searched by varions fnture cosmic-ray measnrements snch 
as CTA |36] and GAMMA-400 [37] . 
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Appendix A: The potential and tree-level annihilation matrices for EWDMs 


Considering the covariant derivative + ig 2 A^T^ for each gange boson A = 

7 , the potential matrix in Eq. ([^ and the tree-level scattering matrix P^ take the 
following general forms: 


c, (r) = 2 5m, 6,, - [T^] 


2 e 


—mj\r 


with Sm, = m^. — ; 


Tial 


Y^ab __ 

2(1 -F 5AB)m 


2 

DM 




_ XA+XB \ 

where f{xA,XB) = ^ 


_ xa+xb '^‘^ 


Xa + Xb ^ (xa - XbY ... 

1-;-h ---- with Xa 
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(Al) 

(A2) 

^DM 


Here the normalization factor W is 1 or \/2 for the Dirac (charged) or Majorana (nentral) 
two-body state, respectively. In the following, we present the explicit matrix elements for 
the three fermionic EWDM candidates in the lowest mnltiplets. 


• Donblet (Higgsino-like) EWDM 

The potential matrix V and the normalized tree-level scattering matrix P^^ normalized 
by (7ra2/ml)^)f(xA,XB) in the basis of the three states XiXij X 2 X 2 ) given 
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as follows: 
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(A3) 


^WW — TW 

lo 


r^z = 
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72 1 1 
72 1 1 

8(1 - 2 sI,Y 272(1 - 2slrf 272(1 - 2s7)2^ 

272 ( 1 - 25^)2 1 1 

1 1 


/ 
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\^ 272(1 - 25^)2 
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0 0 0 
yO 0 Oy 


SwP — 2s 7)^ 
2c^ 


r — 
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'^l 0 0^ 
0 0 0 
yO 0 Oy 


(A4) 


• Majorana triplet (Wino-like) EWDM 

The potential matrix V and the tree-level scattering matrix Tab normalized by 
{7Ta2/m‘f)jPjf{xA,XB) in the basis of the two states {x~^X ~nre given as follows: 


V 


2 5m_i 




\/2o2> 


c^«2e 


— mu/r 



(AS) 


^ww — 


r7z — 


1 

2 


1 1 72 ' 
(^2 2^ 

2s^c^ 0 
0 0 



(A6) 


Let us note that there are a few discrepancies in factors of the scattering matrices 
(A4) and (A6) compared with the previous results in Ref. |3]. 
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• Hypercharged triplet EWDM 


The potential matrix V and the tree-level scattering matrix T^^ normalized by 
(7ra^/mlj^)f(xA,XB) in the basis of the four states (x++x"", X+X", X?X?, X 2 X 2 ) are 
given as follows: 
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